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Abstract 

We develop a formalism for the resummation of target mass corrections in off-forward two- 
photon amplitudes given by a chronological product of electromagnetic currents, arising in e.g. 
deeply virtual Compton scattering. The method is based on a relation of composite operators 
with a definite twist to harmonic tensors, which form an irreducible representation of the Lorentz 
group. We give an application of the framework for the matrix elements of twist-two operators. 
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1 Hadron mass corrections. 



The leading twist approximation to hard processes in QCD is affected by a number of multi- 
phcative and additive corrections. To the first class one obviously attributes radiative corrections 
in the strong coupling, while the latter encodes higher twist contributions, which provide power 
suppressed corrections. In the case when a given large scale Q, that controls the factorization of 
a reaction, becomes rather low, one has to take care of power suppressed effects since they can 
significantly modify the scaling behaviour and the magnitude of the leading twist prediction for 
the corresponding cross section. Power corrections can be divided into two classes according to 
their origin: dynamical and kinematical ones. The first one comprises of multiparton correlations 
inside hadrons and gives rise to new non-perturbative functions. The second one arises from 
a separation of composite operators into parts that have definite transformation properties w.r.t. 
the Lorentz group and thus possess a well defined geometrical twist, i.e. dimension minus spin. 
This decomposition provides the so-called Wandzura-Wilczek type contributions arising from 
components having definite symmetry properties as well as target mass corrections |0] stemming 
from the subtraction of trace terms in the afore mentioned operators. 

In the present paper we develop a formalism for the resummation of the target mass corrections 
for off-forward two-photon processes whose amplitudes are given by a chronological product of 
electromagnetic currents. This method is indispensable for a study of hadron mass effects in the 
cross section of e.g. deeply virtual Compton scattering [|, H, 0, — the processes sensitive to the 
generalized parton distributions (GPDs). The present day facilities can offer the of the order 
of a few (2-4) GeV^ (for the reactions in question) and one has to account for power suppressed 
contributions to the twist-two and -three observables arising from the ratio of the nucleon mass 
to the hard momentum transfer. Unfortunately, the formalism for resummation of target mass 
corrections invented by Georgi and Politzer |Q for the deep inelastic scattering (see also [0 for a 
first discussion of the topic), and accepted in all consequent generalizations and applications to 
forward electroweak scattering amplitudes, is not directly applicable to off-forward processes in 
the context of GPDs. The reason for that is the appearance of towers of new Lorentz structures in 
the matrix elements of local operators, arisen due to the non-zero t-channel momentum transfer, 
which cannot be handled in the fashion proposed in ||^. Note that a discussion of mass effects in 
exclusive processes with a much simpler kinematics can be found in Ref. 

In the present context, instead of dealing with single variable functions, e.g. GPDs, we use 
spectral functions, the so-called double distributions (DDs) ^, which depend on two variables. 
This allows to resum the mass corrections in a straightforward manner and to derive all-order 
results in a compact form. After this is accomplished one can in principle use an inverse integral 
representation to express DDs by means of GPDs, however, this transformation requires the 
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knowledge of the support of GPDs including an unphysical region. Therefore, this transformation 
can be hardly used for a successful phenomenological application. 

The consequent presentation runs as follows. In the next section we give a short introduction to 
a procedure of the trace subtraction by means of harmonic polynomials. In Section |^ we resum the 
target mass effects resulting from the trace subtraction in matrix elements of twist-two operators 
for zero and spin-i targets. Then it is shown that the expansion of these results in powers of 
[M"^ jc^y can be cast in a conventional form of GPDs and be used for numerical predictions of 
target mass effects in physical cross sections. Finally, we conclude. 

2 Harmonic polynomials and twist decomposition. 

The goal of our study is an evaluation of the off-forward Compton amplitude keeping target 
mass effects stemming from the twist-two contributions in the light-cone expansion of the off- 
forward matrix element of the chronological product T),,^(|,— |) = f )| |Pi) of 
electromagnetic currents = ip'jfi'ip- Since these are entirely kinematical effects, we consider 
dynamical higher twist effects given by multiparticle operators as independent and, therefore, we 
safely neglect them. Note, however, that the trace subtraction in twist-three contributions in the 
Wandzura-Wilczek approximation will generate kinematical mass effects accompanied by leading 
twist operators. We will comment on this issue in the concluding section. In the approximation, 
we are considering, the Fourier transformed two-photon amplitude reads 

rri f 74 ix-Q fx x\ ^ f rl^ (J 

i^, = J axe 1^, —2) =-2]^^ + z0]2 

X (P2I cos {x ■ q) 5^.;p/0,(f , -f ) + sin (x ■ q) e^.p/Op(f , -f ) |Pi), (1) 

where S^^.p^ = g ^pQva + 9 ^laQvp — 9 ^ivQ pa and for the totally antisymmetric tensor we use the normal- 
ization e^^^^ = 1. The vector g is a semi-sum of incoming gij, and outgoing g2/x photon momenta, 
1 = ^(^i + 12)- The non-local operators in Eq. are ^Op{x, —x) = ilj{—x)[—x,x]rpip{x) with a 
Dirac matrix Fp which equals 7p or 7p75 for vector {V) or axial-vector (A) sectors, respectively. 
The path ordered exponential [—x, x] will be dropped everywhere in the following formulae. 

We are aiming now in a decomposition of Op into traceless operators TZ^p and TZ^ with definite 
geometrical twist-two and -three, respectively: 

Opix, -x) = nlix, -x) + n%x, -x) + 7^;(x, -x) . (2) 

To this end, we expand first the non-local operator Op(x, —x) in a Taylor series 

Op{x, -x) = 22 — • ■■XpPpm-H^ '^i*^ ^pm-H = i'TpiVf,^ ■ --iVp^ i>, (3) 
j=o J- 



2 



where = Vfj. — T>^ with the covariant derivative = — igB^, and then we perform the 

An effective way of trace 



X 



X O 



twist decomposition for the local operators Op-j 
subtraction is to solve the sufficient condition d^W (x^, d"^) Tj{x) = |T0|, |lT|], where Tj results 
from the contraction of a rank-j tensor with the product of four- vectors . . . Xfj_. , and is a 
subtraction operator. The solution is given by the so-called harmonic polynomials [§] 

H^(x',d'' 



r (j - + 1) 



-1 id^ 



(4) 



For operators having a free vector index one has an extra condition so that in total the trace- 
lessness implies d'^TVp-^ = and dpTV^.j = 0, for i = 2,3, where the twist-two operators are com- 
pletely symmetric while the twist-three operators are antisymmetric in a pair of their Lorentz 
indices. Employing the equation of motion for the massless Dirac field, one can show that 
TZ^ is proportional to total derivative, V, squared, apart from multiparticle operators, 7?,^ oc 
— V^T?.^ + multiparticle operators. When it is sandwiched between hadronic states, it gener- 
ates contributions proportional either to the t-channel momentum transfer squared or given by 
multi-parton correlations functions. Both of them will be neglected in the following. Thus, our 
considerations will be limited by the condition ^ A^, with = = P2 and A^, being the 
hadron mass and t-channel momentum transfer, respectively. 

A straightforward algebra leads to the desired decomposition in operators with definite geo- 
metrical twist: 



P\3 



3 + 1 



TV 



where the traceless twist-two and -three operators read (cf. 

7^L■ = -^dpW+^ {x\d'') x,o„,, 



P\0 



V? 



— \gp„x ■ d - x^dp] 
2j 



1 



9o 



J + 1 



J>0, 

J>1, 



(5) 

(6) 
(7) 



and the remainder is 



TV 



J + 1 



XpX ■ d 



X d. 



x\d' 



It is obvious that the twist-two operators (10) are traceless. Using d^H^ {x^, d"^) O^^- = and Euler 



theorem x ■ d dpH^ (x^, (9^) O^-j 
operators (0) are traceless too. 



(j — l)dpH^ (x^, (9^) Op-j, it is an easy task to check that the 



3 



3 Resummation within DD formalism. 



In this section we will give a general framework for the resummation of kinematical mass cor- 
rections. For simplicity we consider here only the twist-two sector. The major steps consist of: 
a parametrization of the off-forward matrix elements of symmetric local operators by means of 
moments of DDs, the trace subtraction with the harmonic projectors , then a Fourier transfor- 
mation to the momentum space and, finally, a resummation of infinite series. 

In order to demonstrate the main features of our formalism, let us consider first matrix elements 
of twist-two operators sandwiched between states of spinless target (P2|'^p;j|-Pi)- Here only the 
matrix elements of the parity even sector are relevant and read 

(^2rO^^°J...,jPi) = PiA. ■ ■ ■ Ph}Bhw + ■ ■ ■ + ^{A,^ ■ ■ ■ A,,}i?,+i,o + . . . . (8) 

They are build from the two vectors P = Pi + P2 and A = P2 — -Pi, {■ ■ ■} denotes symmetrization 
of corresponding indices and the ellipsis stands for possible terms containing the metric tensor. 
Since we act on this expression with we can neglect the latter. After projection with harmonic 
polynomials via Eq. one gets 



k=0 



Xni ■ ■ ■ Xu. 

Pfll...flj 

with the operation S standing for symmetrization and trace subtraction. The first line gives a 
compact expression for subtracted operator which is extremely convenient for consequent consid- 
erations and is the basis of our formalism. 

On the other hand we can deduce the twist-two operator from Eq. by a contraction of all 
Lorentz indices with a light-like vector n so that the trace terms vanish identically. This tells us 
that the coefficients in front of the Lorentz tensors are related to the reduced matrix 

elements of light-cone operators and can be represented in terms of moments of conventional 
leading twist GPDs or DDs via 

1 d'' 



B 



k\ drj^ 



r]=0 



f^dxx^-'B{x,v) = [I) Jjydzy^-h'fiy,z) (10) 



where 0<k<j,l<j and the integration domain is i7 = { — 1 < y < 1, —1 + \y\ < z < 1 — \y\}. 

Consequently, the parametrizationQ of the matrix element for the non-local symmetric operator 

^This is not the definition given in Refs. ^ which suffered from inconsistencies due to omitted Lorentz struc- 
tures as noted in . Our solution of this problem differs also from where a two-component parametrization 
has been suggested. 
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and the light-ray operators coincide up to terms proportional to Xp or x'^: 

(P2ror(a:,-a:)|Pi) = j ^dy dz f{y, z)Vp e-^^'"" + . . . , (11) 

where we introduced a new shorthand notation = yP^ + -zA^. 

The resummation of matrix elements of traceless local operators can be performed in the 
formalism of DD /(y, z). By means of Eq. (jTUp it is easy to see that the finite sum w.r.t. k in Eq. 
(H) provides 

(AlX^l^i) = X ^^^^ ^) JTT^"^'^' ■ ■ ^^^^ 

The evaluation of projection of harmonic polynomials results into conventional Chebyshev poly- 
nomials related to the irreducible representation of the orthogonal group 5*0(4) 



This result when plugged into Eq. (12) gives our final expression for the twist-two traceless local 
operators. We will use it after we perform the Fourier transformation in Eq. ([l|). We accomplish 
it making use of the following general formula 

where on the r.h.s. we dropped all terms involving the metric tensor g^i^,^-, since they vanish when 
contracted with a traceless tensor. Thus, we get 

= v^^nj + i)j,n,,.'f. . . ^{AIX;....«I'P.). (15) 

with projector 

2 

nf,u = g^,u - -^Qfilu, n^pUpy = Qf^y. (16) 
Then the Fourier transformation of the twist-two contribution to the Compton amplitude gives 

/'''^'=''"[^(-p^i'''^Uf--f)i^i> 

= 2"' htiP^\"n.,W - ln„^tj{PrnlJpA (17) 



Here we used a new convention 



3=0 



_ %1_ ^-7?2 /10\ 
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and exploited an identity 



To evaluate the sums in ([T7|) we substitute ([T3| ) in (0) and replace by q^/q^- This gives the 
local traceless matrix elements (P2|^^p jl-Pi) in terms of Chebyshev polynomials. The summation 
in Eq. ([T^ can be now done making use of the generating function for the Chebyshev polynomials 
Gu{ci, h) = Z^jlo ^^Uji^) = (1 ~ 2a6 + a^)~^. To get a required j-dependent coefficient in the series 
we do integration(s) of both sides with appropriate weight. The sums required in the calculation 
of (0) read 



Y- rf/,+i(6) = / da' 

+ ^ ^ Jo Jo {a") 



iGuia",b) -l-2a%) 



Here a = \\jV'^ /q^ and b = q ■ V j \/q^V'^. This immediately leads to 

YiP^l^'nUP,) = q'nprTT- / dydzf{y,z) (20) 
j—Q Oq^ JQ 

x|^--™infl+^--^™V^\-^™lnfl+^ + -^™n, 



2Vl \^ 2H ; V 2S 

where we introduced the conventions 



M' = — (22) 



q-V (q-V) 
Inserting our findings (pO]) and (|2l|) into Eq. ([ITD the next step is a mere differentiation w.r.t. 
g^, which is done by the formula 



(23) 



for a test function r. Finally, we contract the resulting equation with the tensor S^y-p^, which leads 
to the hadronic tensor 

Tl. = J^dydzfiy, z) {(gV " ^^u) -^i + {v, + |) {v, + |) -^2} , (24) 
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with mass-dependent coefficient functions 



' {1 + M^)[AS{1 + S)-M^]^ 4(1 + ^2)5/2 ''[i + ^/T^M2 + 2s) ^ ' 

In order to get the right iO-prescription, so as to pick up a correct sheet of the Riemann surface 
for the logarithm, we notice that in ([ISD we have to restore the suppressed Feynman prescription 
as follows g2 — > g2 _|_ 

Let us comment on our result. As we observe, the leading order massless (generalized) Callan- 
Gross-type relation (see [0) JF2 = S'^J^i is violated by target mass corrections. We checked 
that our result coincides in the forward limit with the well-known ones j^. Obviously, we have 
QfiT^i, = and therefore current conservation is fulfilled in the forward case. Note that this is only 
the case in the sum of the direct and the crossed amplitudes. However, the gauge invariance in the 
off-forward kinematics is violated, i.e. q2^iT^u ^^id, as it was previously studied in the case of 
massless amplitudes, it is restored once higher twist corrections are accounted for [0, |T^, |TB], . 



This must also persist to the case of traceless operators which are not projected on the light cone. 

Now let us generalize the formalism to a spin-| target. The matrix elements {P2\TZp.j\Pi) are 
build now from the vectors P, A and a set of (independent) Dirac bilinears, characterizing the 
spin content of the target. For our purposes it is convenient to use the following structures 

{h„h,)=UiP2h,il,^,)UiP,), (6,6) =f/(P2) (1,75) f/(Pi), (26) 

for (vector, axial) sectors, respectively. Thus, we parametrize the matrix element of e.g. vector 
operator according to 



(Ar0^^"J...,jPi) = h{,P,,---P,^}A,^^,,+, + ... + h{,A,,...A,^^A,^,,, (27) 

+ ^ {P{pP)^i ■ ■ ■ P^yBj+ij+i + ■■■ + A{pA^^ ■ ■ ■ A^^.}Pj+i,o} + . . . , 

where again terms proportional to the metric tensor are not needed. The non-local twist-two 
vector operator reads 



(2J 



,=oj!(j + l)' V 4 ) '^'\V^J' 

Here we have generated the factor of x in /i ■ x by a differentiation w.r.t. V. The same we did with 
the b form factor for the purpose of a uniform representation, as well as to have a cross check on 
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our previous calculations. Obviously, this step is not required at all and we can readily borrow 
the results already found for a (pseudo) scalar target. The same equation holds for axial operator 
with a trivial dressing of symbols with tildes. 

The reduced matrix elements -Bj+i j+i-fe are represented in terms of GPDs or DDs as in Eq. 
(PUD, — replace f{y, z) by z). For the reduced matrix elements A^+i we write 

1 



A, 



ri=0 



'^dxx^-'A{x,r]) = [\^) j^dydzy^~^-^z^fA{y^z), (29) 



— / dxx^~^ {H{x,ri) + E{x,ri)} = Q. 

UTl^ J — 1 



where j > 1 and the index k now varies in the intervals < A; < j — 1. There is an immediate 
consequence of the fact that j-th moment of the function A{x, rj) is a polynomial of order t]^~^ 
only. Namely, translating it to the conventional Ji's distributions, H = A + B, E = —B, we 
immediately find that 

d_ /■! 

As a particular example we recall that for the j = 2 moment, which corresponds to Ji's sum rule, 
the 77-dependence drops off as it was noted before For the axial channel the GPDs A and B 
are identical to the conventional H and E and this implies that the j-th moment of H (E) has 
the expansion to order j — 1 (j) in r^. 

The Fourier transform and the resummation are done in the same way as before and result in 

= -\jjydzl (fA{y, z)h-d'^ + My, ^) ^ ^ " dA te,,,^q,V, (30) 



+ (f^[y^z)h.d'P + fB{y,z)^V-d'^' 



with massive coefficient functions 



(1 + M'f ~ 2{1+M^)~ 2{1 + M'f + ^ ^ ' 
' ~ 4(l+A^2)5/2 2{1 + M^f ~ 2(l+A<2)5/2 

' 2(l+A^2)i/2 2{1 + M^) 2(1+^2)3/2 I- J 

where 



, 1-VTTM^+2S , , 1 + VT+M^ + 2S 
^± = m ± m 
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Here Li2 is the Euler dilogarithm Li2(x) = — Jo ^ ln(l — y). 

As we mentioned before, the resummed mass corrections can not be expressed in terms of 
GPDs (at least we did not succeed in doing this). However, every term in the mass expansion, 
M^/g^, can be converted into the conventional GPDs representation. Let us demonstrate it for 
the expanded Compton form factor Fi for a (pseudo) scalar target given in Eq. (^). To the first 
non-trivial order, i.e. to O (M^/g^) accuracy, the latter reads in DD form 



F, = dy dz f{y, z) | [cf^ {s-^ - zO) + Cf^ ( 



(o)r_"-i 



ic^'^ {^-' - ^o) + c« {-s-^ - .0) 



9' 

where in our approximation, ^ A^, we set M.^ = —A^y'^S'^ and introduced the coefficient 
functions Ci°\s-^) = - {1 + S)'^ and C[^\s-^) = (1 + + 2 In (^j^) . To cast the DD into 
the usual GPD representations, which are related by the equation x dy dz f{y, z)6{x — y — rjz) = 
B{x,T]), we need the following general result 

/ dydzl'l] J^(s-')fiy,z)= fdxri^ f dx' vt\x.x')[ (^) ^Bix'^r,), 

(32) 

where we used a shorthand notation for the differential operator d{x, rf) = x-^ — rj-^ and the 
kernels 

Vi{x,x) = — ——0^^{x,x-x), With 0^^{x,y) = . (33) 

[n — 1)1 X — y 



Substitution of these expressions back into Eq. ( P^ results into 

Fi = J^dx J^dx' I (c[°'> (l - io) + (-| - io) ^ 6{x - x') (34) 

Cf) (l _ ^o) + ) (-f - zO) ) \//') (x, x') {d{x', f]) - l) d{x', 7])]b{x', r/). 



which has now the desired form. Completely analogous manipulations produce the GPD repre- 
sentation for all other form factors. 



4 Conclusions. 

In the present note, we have presented a formalism for the resummation of target mass corrections 
stemming from the twist-two approximation to the off-forward Compton scattering amplitude. 
The machinery is based on the group-theoretical content of the procedure of trace subtraction 
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from local operators with definite symmetry property - the operators with a well defined twist are 
given by a harmonic projection of the latter. This allows to resum target mass corrections to all 
orders in a compact form in terms of DDs. However, once the amplitudes are expanded in M^/g^ 
series, they can be cast in the form involving GPDs. 

An important feature of the result is that the mass corrections enter with additional powers 
of the scaling variable ^. Thus, we expect that they are suppressed in the small ^ (a^Bj) region. 
Since there still might be a numerical enhancement due to coefficient functions, a numerical model 
dependent study of the target mass effects to the deeply virtual Compton scattering process for 
the kinematics of HERMES and Jefferson Lab experiments is necessary. 

In the twist-two approximation the current conservation is violated and must be restored by 
taking into account contributions with higher geometrical twists. To make the former manifest, 
one has to decompose those higher twist operators in a linear independent set by means of the 



QCD equation of motion |T^, |16|, |T^. This procedure gives us a part which is expressed in terms 
of total derivatives acting on operators with a lower twist. Of course, the latter when combined 
together with the leading prediction must render a manifestly gauge invariant result. However, 
the trace subtraction and the application of the QCD equation of motion do not commute and 
one can easily fail to reproduce an amplitude which respects current conservation. This problem 
goes beyond the scope of the present study and requires further investigations. 



We would like to thank A.V. Radyushkin for a clarifying discussion of Ref. |17 . 
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